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We study the surface critical behavior of branching-annihilating random walks with an even 
number of offspring (BARW) and directed percolation (DP) using a variety of theoretical techniques. 
Above the upper critical dimensions d c , with d c — 4 (DP) and d c = 2 (BARW), we use mean field 
theory to analyze the surface phase diagrams using the standard classification into ordinary, special, 
surface, and extraordinary transitions. For the case of BARW, at or below the upper critical 
dimension d < d c , we use field theoretic methods to study the effects of fluctuations. As in the bulk, 
the field theory suffers from technical difficulties associated with the presence of a second critical 
dimension. However, we are still able to analyze the phase diagrams for BARW in d = f,2, which 
turn out to be very different from their mean field analog. Furthermore, for the case of BARW only 
(and not for DP), we find two independent surface f3i exponents in d — 1, arising from two distinct 
definitions of the order parameter. Using an exact duality transformation on a lattice BARW model 
in d = f , we uncover a relationship between these two surface /3j exponents at the ordinary and 
special transitions. Many of our predictions are supported using Monte-Carlo simulations of two 
different models belonging to the BARW universality class. 

PACS numbers: 05.40.-a, 64.60.Ak, 64.60. Ht 



I. INTRODUCTION 

The study of surface critical behavior in equilibrium 
statistical mechanics has established the importance of 
boundaries in critical systems and their impact on scal- 
ing and universality yj. Quantities measured close to 
the surface can scale differently than in the bulk and can 
possess distinct critical surface exponents. Depending 
on the boundary conditions, various surface universality 
classes are possible each with different values for the sur- 
face exponents. In this paper we will be interested in 
the surface critical behavior of certain dynamic systems, 
which possess a non-equilibrium phase transition from 
an active into an absorbing state from which the system 
cannot escape. 

The most prominent example of a system with an ab- 
sorbing state is directed percolation (DP). It describes the 
directed growth of a cluster governed by a growth proba- 
bility p of its fundamental constituents. For probabilities 
below a critical value, p < p c , the cluster will die after 
a finite time, which means that the system gets trapped 
in the vacuum - the unique empty state. On the other 
hand, for high enough growth probabilities p > p c , there 
is a finite probability that the cluster will always remain 
active. Exactly at p = p c , there is a critical phase transi- 
tion from the active into the absorbing state g| . A whole 
range of other systems possessing a phase transition from 
a non-trivial active phase into a unique absorbing state 
fall into this universality class. Some examples include 
epidemics, chemical reactions, catalysis, and the contact 
process (see || and references therein) . 

During the last few years, however, studies have also 



been carried out for systems with absorbing states which 
do not belong to the DP class. For instance, the model of 
branching- annihilating random walks with an even num- 
ber of offspring (BARW) exhibits quite different behav- 
ior and defines a separate universality class. Other 
models in this class (at least in d = 1) include certain 
probabilistic cellular automata , monomer-dimer mod- 
els [p|-|lO[], non-equilibrium kinetic Ising models |ll[], and 
generalized DP with two absorbing states (DP2) [jL2j. 
These models escape from the DP universality class by 
possessing an extra conservation law or symmetry: for 
the BARW model, a "parity" conservation of the total 
number of particles modulo 2; for the other models, an 
underlying symmetry between their absorbing states. 

In the present paper we study the impact of surfaces 
on the critical behavior of the DP and BARW models. 
Previous work has concentrated on surface effects in DP 
using field theoretic methods jl^Jl^], Monte-Carlo sim- 
ulations in d — 1,2 p^dj^l, the Density Matrix Renor- 
malization Group in d = 1 [ p"7[ , and series expansions 
in d = 1 Relations between surface DP and lo- 

cal persistence probabilities have been explored in [B0[ . 
Work has also been performed on active, but slanted, 
walls in DP, which give rise to a "curtain" of activity 
whose width is given by an angle-dependent correction to 
bulk DP ^l). Critical surface effects for a model in the 
BARW universality class were first studied in . In this 
paper we will build on this earlier work by presenting a 
unified picture of surface critical behavior of both BARW 
and DP. After summarizing the main details of DP and 
BARW in the bulk, we present a comprehensive analy- 
sis of the surface critical behavior of both models using 
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mean field theory. This involves the usual classification 
into ordinary, special, surface, and extraordinary transi- 
tions. However, below their respective upper critical di- 
mensions, fluctuation effects become important in both 
models, and this leads to the breakdown of mean field 
theory. In order to understand this fluctuation regime 
we employ a variety of theoretical techniques. First, we 
construct a phenomenological scaling theory which is able 
to describe the various surface universality classes. This 
scaling theory for BARW can then, to some extent, be 
justified using field theoretic techniques (a field theory for 
surface DP has already been presented in Jl3| ). However, 
the BARW field theory suffers from technical problems 
associated with the presence of a second critical dimen- 
sion, which means that the interesting d — 1 regime can- 
not be accessed satisfactorily. Nevertheless some results 
can be derived field theoretically which we put together 
to draw up a 1 + 1 dimensional surface phase diagram. 
This phase diagram displays many differences from its 
mean field analog. In addition, using exact techniques 
involving a mapping to a quantum spin Hamiltonian, we 
have been able to establish an exact duality transforma- 
tion for a lattice BARW model in d = 1. We find that 
this links together two of the boundary phase transitions 
in d = 1 in a non-trivial way (as was suggested in p2||). 
We have also performed extensive Monte-Carlo simula- 
tions for BARW and DP2 which support many of our 
theoretical conclusions. 

The paper is organized as follows: in Section |H| wc 
briefly introduce the bulk DP and BARW models. Then 
in Section HI we consider the surface behavior of DP, 
where we present an extensive mean field analysis. We 
also summarize detai ls of t he fl uctua tion regime for d < 
d c = 4. In Sections [VA and IV B, we give the phase 
diagram for surface BARW in mean field theory and in 
d = 1 respectively. These results can be con tained within 
a scaling theory presented in Section IV C . We then dis- 
cuss in Section IV D how this analysis can be par tially 
justified using field theoretic methods. In Section IV E , 
we give some exact results for d = 1. Our theoretical 
analysis is then backed up using computer simulations of 
the lattice models introduced in Section M. Details of 
these si mula tions are presented in Section lVE Finally in 
Section VII we round off with some conclusions. 



II. BULK DP AND BARW 



We begin by briefly reviewing the definitions of the 
DP and BARW models. The update rules for bond DP 
in d + 1 dimensions on a tilted square lattice are easily 
defined: for each site at time t, form bonds with prob- 
ability p to the neighboring sites at time t + 1 ||. An 
example of a cluster grown from a single seed according 
to these rules is shown in Figure |l|a. 

For growth probabilities below a certain threshold such 
a process will eventually die out, whereas for higher val- 



ues there is a finite probability of survival, which means 
that the system is in the active state fl23|| . As is well 
known p4]-p6| , various reaction-diffusion models also fall 
into the DP universality class. The simplest of these 
is defined by the following reaction scheme for a single 
species of diffusing particles: 



A^A- 
A + A- 
A -> 



A with rate a 
A with rate A 
with rate /x, 



(1) 



where, in the corresponding lattice model, we allow for 
multiple (bosonic) occupancy of any given site. 



(b) 




FIG. 1. DP clusters grown from a single seed (a) in the 
bulk and (b) next to a wall. 

The second system which we will analyze in detail 
is the BARW model This is defined again by a 

(bosonic) particle model, with the following reaction pro- 
cesses 



A — > (to + I) A with rate a m 
A + A -> with rate A. 



(2) 



For m odd, the above model is known to belong to the 
DP universality class, however for to even we have a new 
universality class. Unless otherwise specified when we re- 
fer to the BARW model we will be referring to the even 
to case. 

The growth of both BARW and DP clusters in the bulk 
close to criticality can be summarized by a set of inde- 
pendent exponents. A natural choice is to consider u± 
and i^i which describe the divergence of the correlation 
lengths in space, £j_ ~ |A| _l/i , and time £|| ~ |A| _I/ ii. 
Here the parameter A describes the deviation from the 
critical point (in mean field theory A = /i — a for DP, but 
A = — ma m for BARW). We also need the order param- 
eter exponent (5, which can be defined in two a priori 
different ways: it is either governed by the percolation 
probability (the probability that a cluster grown from a 
finite seed never dies), 



P(t - co, A) ~ |A| Al 



A < 0, 



(3) 
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or by the coarse-grained density of active sites in the 
steady state, 



»(A)~|A|A»— , A<0. 



(4) 




(b) 



o„, 



When A < the system is said to be in an active state, 
whereas for A = the system is critical (with an alge- 
braically decaying density), and for A > (if applicable) 
the system is inactive (with an exponentially decaying 
density) |27| ]. For the case of DP, it is known that (3 
is unique: /3 see d = /3dens m any dimension, both above 
and below the upper critical dimension d c — 4. This fol- 
lows from field theoretic considerations [^3|,^4) and has 
been verified by extensive numerical work. The relation 
also holds for BARW in 1 + 1 dimensions, a result first 
suggested by numerics and now backed up by an exact 
duality mapping ps|] . However, this exponent equality 
is certainly not always true: if we consider the BARW 
mean- field regime valid for spatial dimensions d > d c = 2, 
then the system is in a critical inactive state only for a 
zero branching rate, where the density decays away as a 
power law. However, any non-zero branching rate results 
in an active state, with a non-zero steady state density 
(see Figure ||a) ||. This density (||) approaches zero 
continuously (as the branching rate is reduced towards 
zero) with the mean-field exponent Pdens = 1. Neverthe- 
less, for d > 2, the survival probability (g) of a particle 
cluster will be finite for any value of the branching rate, 
implying that /3 SO od — in mean-field theory. This re- 
sult follows from the non-recurrence of random walks in 
d > 2. 

Field theoretically, DP is believed to be satisfacto- 
rily understood — the appropriate field theory (sometimes 
called Reggeon Field Theory) J24|,^6| is well under con- 
trol and the exponents have been computed to two loop 
order in an e = 4 — d expansion |p5| . However, for the 
case of BARW, a description of the 1 + 1 dimensional 
case poses considerable difficulties for the field theory 
H . These stem from the presence of two critical dimen- 
sions: d c = 2 (above which mean-field theory applies) 
and d' c w 4/3. For d > d' c the behavior of Figure |^a 
holds, i.e. an active state results for any non-zero value 
of the branching o~ m , whereas for d < d' c the system is 
only active for a m > cr.m.criticai, as shown in Figure ||b 
H . This means that the physical spatial dimension d = 1 
cannot be accessed using epsilon expansions down from 
the upper critical dimension d c = 2. Furthermore, for 
the <7 m < o~ rn . critical region, the system is not inactive (in 
the sense of an exponentially decaying density). Instead 
this entire phase is controlled by the annihilation fixed 
point of the A + A — > process, where the density decays 
away as a power law. Hence this phase should rather be 
considered as still being critical. 



FIG. 2. Schematic bulk behavior for BARW of the density 
n as a function of the branching rate a m for (a) d > 2 and 
(b) d=l. 

Despite the problems associated with BARW for d < 
d' c , we can still put forward a general scaling theory for 
DP and BARW, valid both above and below their critical 
dimensions. However, we must retain a possible distinc- 
tion between /3 SC cd and /3dcns- For example, the average 
lifetime (t) of finite clusters can be derived from the scal- 
ing form for the survival probability 

p(f,A) = |A|a~v(t/e„). (5) 

We then find (t) ~ |A|~ T , where r = v y — /3 SC cd- The ap- 
propriate scaling form for the density n(x, f), given that 
the cluster was started at x = 0, t = 0, is 

n(x,t,A) = |A|^+^-/(z/a,*/£||). (6) 

Notice that rotational symmetry about the seeding point 
x = implies that the spatial coordinates enter the scal- 
ing function only as x — |x|, the distance from the seeding 
point. Using the expression (o) we see that the average 
mass of finite clusters, (s) ~ | A|~ 7 , is related to the other 
exponents via the following hyperscaling relation: 

v\\ + dv± = (3 sccd + (3 dcns + 7- (7) 

Note that (0) is consistent with the distinct upper critical 
dimensions for BARW and DP. Using the above mean- 
field values for BARW and i/± — 1/2, i>\\ = 1, and 7=1, 
we verify d c — 2. In contrast, for DP one has the mean- 
field exponents [3 dcns — /3 scc d = 1 and d c = 4. 



III. SURFACE DP 



We now briefly review the surface critical behavior of 
DP and indicate how the above relations and exponents 
are modified in a semi-infinite geometry, where we place 
a wall at x±_ = [x = (xii,:e_i_), with the _!_ and || direc- 
tions being relative to the wall] . An example of such a 
cluster grown close to a wall is shown in Figure |l|b. 
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FIG. 3. Schematic mean field phase diagram for surface 
DP. See text for an explanation of the labeling. 

A schematic phase diagram for surface DP is shown in 
Figured ( see Ml)i where A s is the deviation of the sur- 
face from criticality. In Figure ^, the labeling conforms 
to the standard nomenclature of surface critical phenom- 
ena: O stands for the ordinary transition (bulk critical, 
surface inactive); Sp is for the special Transition (bulk 
and surface both critical); S is for the surface transition 
(surface critical, bulk inactive); and finally E stands for 
the extraordinary transition (surface active, bulk criti- 
cal). 

The bulk exponents are, of course, unchanged by the 
presence of a surface and, furthermore, one can show that 
the correlation length exponents on the surface are also 
the same as in the bulk. Hence, except at the special 
transition, one finds just one extra exponent: the surface 
density exponent /3i,dcns- This is defined from the steady- 
state density at the wall. For example, at the ordinary 
transition, we have 

n(xj_ = 0, A) ~ | Al^ ^- , A < 0. (8) 

On the other hand, at the multicritical special transition, 
one finds two independent surface exponents - a new sur- 
face density exponent, (3f p dcns , and a crossover exponent 
<j)\. In principle one could also allow for a second type 
of surface j3\ exponent, one defined from a survival prob- 
ability for clusters started on the wall. For example, at 
the ordinary transition, we would have 

Pi(t-> oo, A) - |A|^°=«d, A<0. (9) 

However, the surface exponents here show a similar 
pattern to their bulk counterparts and fulfill Pi sccd = 
/3p dcns = Pi , as can be shown by a field-theoretic deriva- 
tion of an appropriate correlation function [fl4f . This kind 
of equality should also hold for the j3\ exponents at the 
special and surface transitions. 

Numerically, the exponents at the ordinary transition 
have been measured very accurately using series expan- 
sions (for d = 1 ) [|18 [ |19| ], and Monte-Carlo simulations 
(for d = 1,2) |l4|-|l6||. However, there has been no numer- 
ical work to date on any of the other possible transitions 
on the boundary. 



A. Mean Field Theory 

Although a considerable amount of work has already 
been performed on surface DP [p^|- p^|p^| , a comprehen- 
sive mean field analysis has been lacking. The purpose 
of this section is to provide such an analysis, and in the 
process we will derive several new results. The equation 
describing mean field DP with a surface is 

d t n = DV 2 n- An-\n 2 , (10) 

with the boundary condition 

Dd x± n\ x±=Q = A s n\ x±=0 . (11) 

Here the variable A = fx — a is the difference between the 
rates for the A — > and A — > A + A processes. Similarly 
we have the surface variable A s , and the bulk quadratic 
term is due to the reaction A + A — > A. Note that a sur- 
face A + A — > A reaction does not have to be included, 
as it is an irrelevant process in the renormalization group 
(RG) sense From the above equation (|l(i|), the bulk 
mean field exponents can easily be computed: u« = 1, 
v±_ = 1/2, and /3 = 1. Furthermore, with the inclusion of 
a boundary, we see that the correlation length exponents 
are unchanged at the wall but the surface f3i exponents 
are altered. If we are interested in the mean field steady 
state, then we can replace ([!(]) with 

Dn" -An- \n 2 = 0, (12) 

where n" = d 2 n/dx 2 L . The appropriate boundary condi- 
tion ( fTl| ) is given by 

Dn' s = A s n s , (13) 

where n s = n\ x± =o, and n' s = dn/dx±\ x±= Q. Multiplying 
( |l2| ) by n' and integrating, we have 

^Dn' 2 - ^An 2 - ^Xn 3 + C = 0, (14) 

where C is a constant of integration. Using the bulk re- 
sults n' = 0, and n — (— A)/A for A < 0, or n — for 
A > 0, we have 

= 1/2 (n. ^ 4 1 ) (§». + W) V2 ^ < 0] (15) 
^ = -[£] 1/a MK + £) Va [A>0](16) 

where we have also used the boundary condition (|l3|). 

• Ordinary Transition. Consider the case where A s > 
and A — * CP. In that case we expect n — |A|/A 3> n s , 
and thus Eq. ((TJ) yields n s cx |A| 3 / 2 , giving $> = 3/2. 

This exponent can also be derived on physical grounds 
as follows (see also p9[). At the ordinary transition the 
density falls to zero not exactly at the wall but would 
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rather reach zero a distance I on the far side of the sur- 
face (if the density were appropriately continued). Hence, 
the density on the boundary can be computed from £n' s . 
Since I is a microscopic distance which remains finite 
even at the bulk critical point, one can compute the scal- 
ing of the surface density simply from n' s = dn/dx±\ x±= Q. 
Thus, from dimensional analysis, we see that =3/2. 

• Special Transition. In this case if A s = we see from 
( |l5| ) that n s scales in the same way as the bulk density 
n, i.e. f3f p = 1 pfj|| . Furthermore, a simple rewriting of 
( |l5| ) and (|l6|) reveals the scaling A s ~ A 1 / 2 , fixing the 
crossover exponent as <f>\ = 1/2. 

• Surface Transition. For this case A s < and A > 0, 
and hence from Eq. © we find n s = (3 /2DX) [A 2 -DA] 
for < DA < A 2 , and n s = for DA > A 2 . Hence 
the line in parameter space where the mean field surface 
transition occurs is given by A 2 = DA cr iticai, and we then 
have Pf — 1. Note that this is the same value as in the 
bulk, a standard feature of the surface transition which 
is believed always to be in the same universality class as 
a d — 1 dimensional bulk transition. As we are dealing 
with mean field theory this will of course yield the same 
exponent for the surface transition as in the bulk. 

• Extraordinary Transition. In this case the surface 
density is of course non-zero both above and below the 
transition. However, if we expand n s in powers of A 
for A s < and A — > + or 0~, we see that these two 
expansions differ at third order plf , i.e. n s has a dis- 
continuity in its third derivative at A — 0. Hence we 
identify 0f = 3. To the best of our knowledge, this tran- 
sition does not seem to have been previously discussed in 
the literature. 

However, as is the case in equilibrium critical phenom- 
ena, we expect the extraordinary transition to be more 
general than the scenario described above. In fact, the 
extraordinary transition is associated with the onset of 
order in the bulk regardless of how the surface is ordered. 
In particular, for arbitrary values of A s , the surface can 
be ordered by applying the equivalent of a surface ex- 
ternal magnetic field. For the BARW process this is 
simply the surface spontaneous particle creation reaction 
— > A. Extending our previous mean field analysis to 
cover this case (sometimes called the normal transition), 
we recover precisely the same results as obtained above, 
with j3f — 3. Hence the important point for the extraor- 
dinary transition (as described in |32| ) is that the surface 
must be active at A = — the means by which this is 
achieved is unimportant. 

Next, we consider the case where the bulk is exactly 
critical, i.e. A is exactly zero, and therefore the corre- 
lation lengths £|| and £j_ diverge. In that the case the 
density in the bulk decays away as dn/dt = —An 2 => n ~ 
t . Hence, for the surface, we need now to include time 
dependence in our analysis, and therefore ( |l2"| ) is replaced 
by 

h = Dn"-Xn 2 , (17) 



where h = d t n and the boundary condition remains as 
given in Eq. (|||). Multiplying both sides of (0) by n' 
and integrating, we obtain 

f°° 11 
J dx ± nn' = -—A 2 s n 2 s --X(n 3 -n 3 s ), (18) 

where we have have used the conditions n' = (in the 
bulk) and the boundary condition (|l3|). 

• Exactly at the Extraordinary Transition. Here, the 
density close to the wall will be in an active steady state, 
and hence nearby n will be close to zero. However, well 
away from the surface we expect to recover bulk behav- 
ior where n' « and n ~ t . Hence, to leading order 
the integral on the LHS of Eq. ( |l8| ) will be zero. There- 
fore, from (p^|), we find a steady state on the surface with 
n s 3A 2 /2DA. Furthermore, we expect that this active 
region will extend into the bulk, with the density decay- 
ing away asymptotically as x^ 2 . However, assuming the 
system is started with initial conditions at t = of con- 
stant density everywhere, then after a time t this region 
will only extend into the bulk as far as x±_ ~ t 1 / 2 , where 
we will find a crossover to the bulk t" 1 density decay. 

• Exactly at the Special Transition. Here, where 
A s = 0, we see that the mean field equation ( |l8| ) is solved 
by n s = n ~ t^ 1 . Hence the surface density scales in the 
same way as in the bulk (see also Q). 

• Exactly at the Ordinary Transition. Once again if we 
start with initial conditions of uniform density at t = 0, 
then at later times a depletion zone will be formed close 
to the surface. This zone will again extend a distance 
of order t 1 / 2 into the bulk. The surface scaling can now 
most simply be derived via dimensional analysis of the 
surface operator d x± n\ x±= o, yielding n s ~ i~ 3 / 2 . 



B. Beyond Mean Field Theory 

We expect that the phase diagram shown in Figure || 
is generally valid for surface DP close to the upper crit- 
ical dimension d c = 4. However, in 1 + 1 dimensions, 
where the surface is just a zero dimensional point, the 
phase diagram may look rather different. For example, 
for an inactive bulk, net particle production is only pos- 
sible at one point. Furthermore, since particles will be 
constantly lost into the bulk, where they will decay away 
exponentially quickly, it will probably not be possible to 
form an active surface state (at least for finite particle 
production rates). If this is the case only the ordinary 
transition will be accessible in d = 1. Furthermore, for 
arbitrary dimension, we note that a system which is sim- 
ply cut off at x±_ — can also only undergo an ordinary 
transition. This is a result of there being the same micro- 
scopic reaction rates on the surface as in the bulk. In low 
dimensions it becomes more and more difficult to induce 
an active state (since the fluctuations become larger), and 
hence if the bulk is adjusted to be at criticality, it follows 
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that the surface (considered independently) would be in- 
active. Therefore, for the case of DP, one will only be 
able to find the ordinary transition (as was certainly the 
case in the simulations of 0]). 

The scaling forms for the survival probability (at the 
ordinary transition) and correlation functions (at the 
special and ordinary transitions) have been discussed in 
|p3| , pr| . For example, at the ordinary transition, the sur- 
vival probability for a cluster started on the wall at t = 
has the scaling form jl4| 

P 1 (t,A) = \Af?cp 1 (t/^). (19) 

Hence the average lifetime of finite clusters at the ordi- 
nary transition, (i) ~ |A| _r i , satisfies t® = Mi — a 
straightforward generalization of the bulk result. Previ- 
ous series expansions in 1 + 1 dimensions have indicated 
a value for t° equal to unity fl8fl , although very recent 
(and even more accurate) series results [|l9| have cast 
some doubt on this conclusion. No theoretical explana- 
tion for why should be equal to unity has emerged. 

The bulk density n(x, t) for a cluster initiated on the 
wall at t = is given by p3| . pT| 

n°(x, t, A) = |A|# [x/Cl, t/$|] . (20) 
n^(z,i,A,A s ) = lAfi^h [x/tL, t/£ p A S /\A\<»] , 

where the surface exponents have been calculated to 
0(e = 4 - d) in p| giving 0? = 3/2 - 7e/48, (3f p = 
1 — e/4, and 4>\ = l/2 — e/16. The first of the expressions 
in ( p0| ) refers to the ordinary transition and the second to 
the special transition. Crudely speaking, the A-prefactor 
in (^o|) comes from (|l9|) for the probability that an infi- 
nite cluster can be grown from the seed, and from (^) for 
the conditional probability that the point (x, t) belongs 
to this cluster. At the ordinary transition, for example, 
it is then straightforward to derive hyperscaling relations 
for the mass of finite clusters which are seeded on the 
wall. This mass scales as (s\) ~ |A| -71 , and as shown 
in (TJ|] one finds 

I/,, + dv±_ = /?°+/? + 7 °. (21) 

For a more detailed discussion of the DP ordinary tran- 
sition, including other types of hyperscaling, we refer to 
Rcf. @. 

IV. SURFACE BARW 

We now turn our attention to the main object of this 
paper, an understanding of the surface critical properties 
of BARW. We will begin by discussing the BARW sur- 
face phase diagram in various dimensions. The basic idea 
is that on the surface we may include not only the usual 
branching and annihilation reactions but potentially also 
a parity symmetry breaking A — ► reaction. Depend- 
ing on whether or not the A — > reaction is actually 



present, we may then expect different surface universal- 
ity classes according to whether the symmetry of the bulk 
is broken or respected at the surface. A similar situation 
in an equilibrium system has recently been analyzed in 
Ref. [B4J. We will find that the competition between the 
parity breaking A — » reaction and the BARW processes 
gives rise to some interesting phase diagrams. 



A. Mean field phase diagram 

The surface phase diagram for the mean field theory 
of BARW (valid for d > d c = 2) is shown in Figure |^. 
Here a m , a mg are the rates for the branching processes 
A — > (m+l)A in the bulk and at the surface, respectively, 
and /i s is the rate for the surface spontaneous annihila- 
tion reaction A — > 0. Otherwise, the labeling is the same 
as that for the DP phase diagram (see Figure ||) . 

The first feature to note is that the bulk is either ac- 
tive (<r m > 0) or critical (a m = 0), but never inactive. 
Hence, unlike DP, there is no possibility of finding a sur- 
face transition, where the surface is critical with the bulk 
inactive. For the case where a m — fi s — 0, we expect that 
for any finite value of the surface branching, the surface 
will become active. This corresponds to the extraordi- 
nary transition with an active surface and critical bulk. 
On the other hand for a m = a mg = and /i s > 0, the 
density at an (isolated) surface would decay away expo- 
nentially quickly due to the A — > reaction. Hence the 
bulk is critical, with the surface inactive; i.e. the ordi- 
nary transition. Consequently with a m = 0, but both fx s 
and a mg non-zero, there should be a line of special tran- 
sitions dividing the extraordinary and ordinary regions. 
This explains the general features of the phase diagram 
in Figure ||. 

At a more quantitative level, the mean field equation 
for BARW is very similar to that for DP: 

d t n = DV 2 n- An -An 2 , (22) 

with the boundary condition 

Dd x± n\ x±=Q = A s n\ x±=0 . (23) 

However the values of the A, A s parameters are now dif- 
ferent: A = —ma m and A s = —mo ms + ji s . The fact 
that A is always non-positive excludes any possibility of a 
surface transition. Otherwise we expect the same mean 
field exponents as in DP for the special and ordinary 
transitions |3(J (except for the /3i )Se ed exponents, see be- 
low). However, the non-positivity of A also leads to an 
ambiguity associated with the definition of (3f dens ; we 
would have to know the behavior of the surface density 
on both sides of the extraordinary transition if we wanted 
to isolate the discontinuity and extract the exponent. 



G 




FIG. 4. Schematic mean field surface phase diagram for 
BARW. See text for an explanation of the labeling. 

We can also consider the mean field behavior of the 
/3i,seed exponents, which is very different from the cor- 
responding behavior in DP. Consider placing two parti- 
cles next to the surface at t = 0. From the recurrence 
properties of random walks we see that, regardless of 
the reaction rates on the surface or in the bulk, there 
is a finite chance that the two particles will never meet 
again. Hence the survival probability is non-zero and 
thus /3i )See d = in mean field theory for the ordinary 
and special transitions. 



B. Phase diagram in I + I dimensions 

Next, we turn our attention to the phase diagram for 
I + I dimensions shown in Figure [5]. Although we will 
make a few remarks below, we will postpone a proper 
justification until we have discussed the appropriate field 
theory in Section IV D. The phase diagram looks quite 



different from its mean field analog due in part to the 
shift of the bulk critical point away from zero branching 
rate, but also due to the absence of any extraordinary 
transition (for finite reaction rates). Physically, this is 
due to the fact that excess particle production (with a fi- 
nite reaction rate) at a zero dimensional surface is simply 
not efficient enough to generate an active state, due to 
leakage into the critical bulk (which for a m < cr mi criticai 
is controlled by the fixed point of the A + A — > reaction 

§)■ 



(a) 



(b) 



O* 



Sp 



o 



FIG. 5. Schematic surface phase diagrams for BARW in 
d — 1 for (a) a m < <r m ,criticai, and (b) <r m = CT m ,critical. See 
text for an explanation of the labeling. 

However, for infinite branching rates and/or if the reac- 
tion — > A is added at the surface, then an extraordinary 
transition should become accessible, although we will not 
consider this case in any further detail. The other main 
features of the phase diagram in 1 + 1 dimensions are 
outlined below. 

• Sp*: For a m < a m , critical and [i s = 0, the system 
is controlled everywhere by the annihilation fixed point. 
In that case one has the special transition, but now in a 
slightly different sense to what we have seen before. In 
this region it is not possible to obtain an active state ei- 
ther on the surface or in the bulk, by small changes in the 
bulk and/or surface branching rates. Hence this "transi- 
tion" is actually entirely controlled by the A+A — > pro- 
cess with the branching playing essentially no role. Thus 
we have marked this "transition" as Sp* in Figure ||a. 
This simpler and analytically tractable case has already 
been extensively analyzed in Ref. |33] . We will po stpone 



further theoretical discussion until Section IV D 4 



• Sp: Next we consider the special transition at a m = 
c m , critical, Ms = 0- This transition borders the bulk active 
phase, and hence will belong to a quite different univer- 
sality class to that described immediately above, and will 
instead be similar to the special transitions discussed in 
earlier sections (although fluctuations will now be very 
important for this 1 + 1 dimensional case). 

• O*: For /i s > and a m < a m> critical, the presence 
of the A — > reaction on the surface gives rise to an 
ordinary "transition" . However, as explained above, the 
branching process again plays essentially no role here. 
Further details of this O* "transition" are provided in 



Section IV D 4 



l-ls 



l-ls 



• O: Finally, at a m = a m , critical, Ms > 0, we expect 
an ordinary transition similar in character to the ordi- 
nary transitions discussed in previous sections (although 
in this 1 + 1 dimensional case the fluctuations are again 
very important). 



C. Scaling Theory 

In this section we construct a scaling theory for the 
survival probabilities and correlation functions at the Sp 
and O transitions. This scaling theory is certainly valid 
for the mean field regime but there are, however, sub- 
tleties involved in its application to the fluctuation dom- 
inated regime for d < 2. The extent of its validity in 
that region will be discussed in detail in the next sec- 
tion. When writing down this scaling theory we must 
also bear in mind the important distinction between the 
/3i,dcns and /3i, scc d exponents. We begin by giving a scal- 
ing form for the survival probability Pi(t, A), where A is 
the deviation from bulk BARW criticality. For example, 
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at the ordinary transition, for a seed placed on the wall 
at x = 0, t = 0, we have 



Pi(t,A) = lAl^cd^/^ ) 



(24) 



It is then straightforward to compute the average lifetime 
of finite clusters, (t) ~ |A|~ ri , where r° = i/y — /?° seod , 
just as in the case of DP. 

Next, we consider the coarse-grained particle density 
m at the point (x, t) for a cluster grown from a seed 
located next to the wall at x = 0, t = 0. At the ordinary 
transition we have 



m(x,t,A) = \A\fi 



,+0i 



s 9i (x/Il, V^ll) • (25) 



As was the case for DP, the A-prefactor in ( p5| ) comes 
from (24) for the probability that an infinite cluster can 
be grown from the seed, and from ([|) for the (conditional) 
probability that the point (x, t) belongs to this cluster. 
The shape of the cluster is governed by the scaling func- 
tion gi and we assume that the density is measured at a 
finite angle away from the wall. If the density is measured 
along the wall, we have instead 



n n (x,t,A) = \Af 



3 5n (x/tt, t/€||) , (26) 



as we pick up a factor lAI^ 1 ^ " 3 rather than lAI^™ 5 from 
the probability that (x, t) at the wall belongs to the clus- 
ter. 

The above correlation functions need only be modi- 
fied slightly to be valid at the special transition. If A s 
(the deviation of the surface from criticality) is a relevant 
parameter, then we must take care to include the extra 
variable As/IA^ 1 in the scaling function. The scaling 
form replacing (|25|) then becomes 

mOM.A) = \AfZ^+^~ gi [ x/u> t/€||; As /|A|^] , 

(27) 

where 4>i is a crossover exponent associated with the mul- 
ticritical special transition. Similarly Eq. ( prjj ) is replaced 
by 

niifot.A) = \AfZ^+^s~ gn [ x /U,t/Z h A s /\A\^]. 

(28) 

Note, however, that there are subtleties concerning the 
special transition in 1 + 1 dimensions which will be dis- 



cussed in Section [V D 1 



At the ordinary transition, for example, we can use 
the above scaling forms to derive some further exponent 
equalities. The average size of finite clusters 



{si) 



"7? 



(29) 



follows from integrating the cluster density (|25|) over 
space and time, where the surface (susceptibility) expo- 
nent 7f is related to the previously defined exponents 
via 



v» + dv±_ = /3? seed + /9dons + 7i°- 



(30) 



Analogously, by integrating the cluster wall density ( |2q ) 
over the (d— l)-dimensional wall and time, we obtain the 
average size of finite clusters on the wall 



IAI-tm. 



where 



I/,, + (d - IK = ^secd + /?i°dc„ S + 7i°i. 



(31) 



(32) 



Note that if the 7 susceptibility exponents obtained from 
(|30| ) and ( |32]) are negative, then they should be replaced 
by zero in (||) and @. 



D. Field Theory 



In order to properly understand the effects of fluctua- 
tions, and to justify some of the scaling forms proposed in 
the last section, we now turn to the development of a field 
theory for surface BARW. We will begin by reviewing the 
field theory for BARW in the bulk, before moving on to 
derive the appropriate surface actions. The bulk BARW 
action, written in terms of the response field ip(x,t) and 
the "density" field tp(x,t), is given by M 



d d x 



dt 



iMx,i)[d t -£V 2 Mx,i) 

-A[l-^(x,i) 2 ]V(x,i) 2 (33) 
+a m [1 - i>(x, t) m ] i>(x, t)^(x, t) 

-ip(x,r) -n ^(x,0) . 



Here the terms on the first line of (|3^) represent diffusion 
of the particles (with continuum diffusion constant D). 
The second line describes the annihilation reaction (with 
continuum rate A), whilst the terms on the third line rep- 
resent the branching process (with continuum rate a m ). 
The final two terms represent, respectively, a contribu- 
tion due to the projection state (see p5fl), and the ini- 
tial condition (an uncorrelated Poisson distribution with 
mean no). In the following we will restrict ourselves to 
the case of even to, since it is known that the odd to case 
belongs to the DP universality class ||. 

The action given in ( |33"|) is a bare action. In order to 
properly include fluctuation effects one must be careful to 
include processes generated by a combination of branch- 
ing and annihilation. In other words in addition to the 
process A — > (to + 1)A, the reactions A — > (to — 1)A, . . . , 
A — > 3A need to be included. These considerations lead 
to the full action 
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Sbulk[V>: V^; t] = I d d x 



dt 



m/2 



+ ^[1 - &x, t) 2l ]^(x, i)V(x, t) (34) 
1=1 

-A[l - $(x, t) 2 ]^(x, i) 2 ] - ^(x, t) - n ^(x, 0)] . 

Notice also that (for even m) the action (|34j) is invariant 
under the "parity" transformation 

<A(x, t) -» -V>(x, t), ^(x, -» -V(x, *)• (35) 

This symmetry corresponds physically to particle con- 
servation modulo 2. The presence of this extra symme- 
try now takes the system away from the DP universality 
class, and into a new class: that of branching-annihilating 
random walks with an even number of offspring. 

Close to the upper critical dimension d c = 2, the renor- 
malization of the above action is quite straightforward 
(here we quote the results from H). Only the branching 
and annihilation rates need be renormalized, and in par- 
ticular there are no diffusion constant or field renormal- 
izations. Furthermore, if we are close to the annihilation 
fixed point, then the RG eigenvalue of the branching pa- 
rameter becomes y„ m = 2 — m(m + l)e/2 + 0(e 2 ), where 
e = 2 — d. Hence we see that the lowest branching pro- 
cess is actually the most relevant. Therefore, close to 2 
dimensions where the branching remains relevant, we ex- 
pect to find an active state even for very small values of 
the branching (in agreement with the mean field phase 
diagram). Furthermore, the fact that there is only one 
eigenvalue to compute perturbatively (as the renormal- 
ization of the annihilation rate can be performed to all 
orders p6|), means that there is only one independent 
exponent. Hence, close to 2 dimensions, the order pa- 
rameter exponent can be shown to be /3dcns = du±. 

However, inspection of the above RG eigenvalue y Gm 
shows that it eventually becomes negative (if the one 
loop result is to be believed). In that case we expect a 
major change in the behavior of the system, since the 
branching process will no longer be relevant at the anni- 
hilation fixed point. The critical transition point is then 
shifted with the active state only being present for values 
of the branching greater than some positive critical value 
(as indicated in Figure ||d). Consequently, we see that 
there is a second critical dimension d' c w 4/3 whose pres- 
ence immediately rules out any possibility of accessing 
the non-trivial behavior expected in d = 1 via perturba- 
tive epsilon expansions down from d = 2. Instead cruder 
techniques (such as the loop expansion in fixed dimen- 
sion) must be employed 

We now turn to the derivation of the surface actions 
appropriate for the cases jj, s = and (i s ^ 0. 



1. (i B — Field Theory 

Starting from an appropriate master equation for the 
system on a lattice, the form of the surface action can 
be derived using standard techniques |35|;[36| . After map- 
ping to the continuum theory, we find the bare action 
+ 5 baxe , with S'buik given by @ and 



cbare cbarc 

° — °bulk 



dt 



(36) 



where ips = ip(x\\,x±_ = 0,i) and ip s — ip(x.n,x_i = 0,t). 
Note that the terms representing the annihilation reac- 
tion A + A — > are irrelevant on the surface close to the 
upper critical dimension.. The classical field equations 
for the above action can be derived by taking the varia- 
tional derivatives of the action 5 baro = 5^ + 5 baro with 
respect to the fields ip and ip. These equations are solved 
by tp = tp s = 1, with satisfying 



(37) 



where A = — ma m , and with the boundary condition 
Dd x± ip\ x±= o = A S -0 S , where A s — —ma ms . These mean 
field resul ts arc in agreement with our analysis in Sec- 
tion IV A. Furthermore, we note that a boundary term 
of the form tp s d x± tp\ x± =o, although marginal from power 
counting arguments, is actually always redundant (even 
in the regime where mean field theory no longer applies). 
This is also the case for the surface action in DP (see 

However, if we are properly to include fluctuation ef- 
fects, we must again take care to include surface terms 
generated by a combination of branching and annihila- 
tion (as in the bulk) . This leads to the full surface action 



Si 



id-l r 



dt 



m/2 

E 

i=i 



o-2i s [i - i>? 



(38) 



Note also that the parity symmetry ( p5|) is preserved for 
the fj, s = model at the wall, as well as in the bulk. 

Power counting on the above action reveals that the 
surface branching rates <Jii a all have naive dimension 
[o~2i s ] ~ where k denotes an inverse length scale. 
However, below 2 dimensions this scaling dimension will 
be renormalized downwards (this can be seen physically 
as a result of processes like A — ► 3A — ► A rendering the 
branching process less efficient). As a result of this renor- 
malization, we expect the lowest generated process (i.e. 
with I = 1 in Eq. (|38|)) will become the most relevant 
(as it was in the bulk). Nevertheless, despite this down- 
ward renormalization, close enough to 2 dimensions, the 
scaling dimension of the most relevant coupling <72 s will 
remain positive, and thus under the RG will flow to co for 
all non-zero starting values. This state of affairs corre- 
sponds to the extraordinary transition, where the surface 
is active while the bulk is critical. On the other hand, at 







bulk criticality and with <72 s =0, we have a multicritical 
special transition point. In this case, after writing down 
and solving the appropriate RG equations (exactly along 
the lines of ^Ju|]), one can der ive th e scaling results for 
the density quoted in Section IV C , where we can see 
that the independent renormalization of <72 s contributes 
to the crossover exponent </>i. Furthermore, since there 
is no field renormalization (either at the surface or in the 
bulk), this implies that the exponent flf dens is just the 

same as in the bulk, i.e. f3 1 dons = /Sdons- However, we 
must again stress that this result is only true close to 
d=2. 

The situation in d = 1 is rather different, partly due 
to the shift of the bulk critical point away from o~ m = 0. 
This means that the d = 1 transition at a m = cr m critica i 
cannot be based on perturbative epsilon expansion calcu- 
lations down from 2 dimensions. However, we can say a 
little more if we first consider the regime a m < o~ m .critical 
in d = 1, where the bulk is controlled solely by the 
A + A — > reaction. In that case we expect the scaling 
dimension of all the a%i s to be negative in d = 1, following 
the downwards trend in the renormalization mentioned 
above. In that case surface branching is then irrelevant in 
d = 1 leading to the Sp* special "transition" . Similarly, 
at the Sp transition at a m = CT m , C riticai, we might again 
expect the <7 2 ;, to be irrelevant. This will be reflected 
in the scaling functions for the density, where (unusu- 
ally) the crossover term A^/IAI^ 1 will now be absent. 
However, the surface exponents here will presumably be 
unrelated to the bulk exponents, since the absence of field 
renormalization mentioned above is not expected to hold 
all the way down to d = 1. 

Hence, if the above scenario is correct, we do not ex- 
pect to see an extraordinary transition in d = 1 for any 
finite value of the surface branching, since the surface 
branching will always be irrelevant. We have confirmed 
this analysis numerically: our simulations have found no 
evidence of an active surface state for <j m < <x m , critical 
even for very high values of the surface branching pa- 
rameter in a fermionic lattice model in 1 + 1 dimensions 



d t $ = DV 2 tp -Aip- 2Xip 2 



(40) 



(see Section VI for further details). 



2. fi s + Field Theory 

In this case the reaction A — » is now possible, but 
only at sites on the wall. In the bosonic field the- 
ory language employed above, we have the bare action 

5 barc = ^bare + ^bare^ where 



£barc = I I dt [a m \\-^]^ a 



(39) 



+fJ,s[*Ps ~ i\lpa 



The symmetry ( p5[ ) is now broken by the surface term 
proportional to /i s , which describes the A — > reac- 
tion. Repeating our derivation of the classical (mean 
field) equation for the ip density field, we find 



where A = — mo~ m , and with the boundary condition 
Dd x± ip\ x±=0 = A s ip s , where A s = -ma ms + fj, s . This 
is in ag reem ent with the mean field analysis given in 
Section IV A . Note that a boundary term of the form 
ip s d x± ip\ x±= o is again always redundant. 

The action ( |3^ ) is a bare action whose terms simply 
represent the reactions A — > (m + 1)A and A — > at 
the surface. Clearly, however, from these two reactions 
we can generate the hierarchy of processes A — > mA, 
A — » (m — 1) A, . . . , A — * 2A. Hence we must replace the 
above bare surface action with 



So = 



dt 



(41) 



.1=1 



The renormalization of the action ([4l| ) is now some- 
what different from the fi s — case. We again expect 
that we need only keep the lowest generated branching 
term on the surface, namely that with I = 1 in (|4l]). 
As before, we expect fluctuations to lower the scaling 
dimension of this coupling from its mean field value (al- 
though actually in d = 2 this suppression will only be 
logarithmic). On the other hand, the efficacy of the 
A — > reaction is certainly not reduced by fluctuations. 
Hence, we expect that A s — > fi s — a% 3 will always run 
to the fixed point at oo, corresponding to the ordinary 
transition. In that case the surface 0° dcns exponent is 
again simply related to the bulk result due to the ab- 
sence of any surface field renormalization. This exponent 
can be computed from the scaling of the surface opera- 
tor Oip / dx ±\ x ±= q , where the x±_ derivative simply brings 
out an extra factor of v± from the scaling function, giv- 
ing /3° dens = /3dens + v^. Using the result ^ d cn S = dv ± 
from Q, we see that /3p dcns = (d + l)u±. Again we stress 
that this result is only true close to 2 dimensions. The 
more interesting transition at o~ m = o~ m , critical in d = 1 is 
not perturbatively accessible in epsilon expansions down 
from d = 2. Nevertheless, we still expect the same gen- 
eral picture to hold with the surface branching always 
being irrelevant, leading to the O* (a m < a m , critical) or 
O (cr m = er m ,criticai) transitions. 



3. Discussion 

From the above analysis we can understand the struc- 
ture of the phase diagram close to 2 dimensions. For 
fi s = this is similar to the mean field picture, with a 
special transition point at a m = a m3 = 0, and with the 
extraordinary transition for a m = 0, o~ ms > 0. On the 
other hand, for fi s > 0, the picture is very different from 
mean field theory, with renormalization effects ensuring 
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that only the ordinary transition is accessible. However, 
actually in d — 2 this might be hard to observe, since 
in that case the surface branching is only marginally less 
relevant. 

One would now like to use the actions ( |36| ) and (|3^) as 
the starting point for a field theoretic investigation of the 
c m = o~ m . critical transitions in d = 1, where one would 
like to identify two independent, non-trivial surface f3\ 
exponents (a feature which is certainly indicated by our 
simulations, see Section VI). Surprisingly, our numeri- 
cal results also indicate that these surface Pi exponents 
"swap" if the /i s = and fi s ^ cases are interchanged 

(i-e- ^dcns = /?i P seed and #Ud = £%J- These inter- 
esting results certainly merit further analysis. Unfortu- 
nately the use of field theoretic techniques here will be 
plagued by precisely the same problems as afflicted the 
bulk calculation, namely the appearance of a second crit- 
ical dimension d' c . Hence one would be forced into using 
uncontrolled techniques (such as the truncated loop ex- 
pansion in fixed dimension) whose values for the bulk ex- 
ponents are known to be in rather poor agreement with 
numerics 0. Furthermore, field and diffusion constant 
renormalizations, which will be of considerable impor- 
tance in d = 1, are not adequately taken into account in 
the truncated loop theory. In fact these renormalizations 
only appear at two loop order. Unfortunately the authors 
of p] were unable to show that a meaningful truncated 
loop theory exists at all at the level of two loops. In 
addition, further technical difficulties exist for d < d' c in- 
volving dangerous irrelevant variables, which have so far 
prevented a derivation of scaling relations at criticality 
even in the bulk. In the light of these problems we have 
not attempted to extend the truncated loop analysis to 
surface BARW. 



4- Results from A + A 



terizations of the Sp* state do not scale in the same way 
close to the wall. 

Since the bulk at the d = 1 O* transition is controlled 
solely by the reaction A + A — > , its properties can 
also be inferred. The relevant surface operator is again 
just d x± ip\x ± =o- Therefore, since distance still scales as 
[a;] ~ i 1 / 2 , we can obtain the required surface scaling 
from simple dimensional analysis. Since the bulk d = 1 
density decays as i -1 / 2 , we see that the surface density 
must decay as t . 

However, we must emphasize again at this point that 
the methods and results mentioned here are only applica- 
ble where both the surface and bulk branching processes 
are unimportant. Unfortunately, therefore, the more in- 
teresting Sp and O transitions in 1 + 1 dimensions remain 
out of reach. 

Hence, given the fundamental difficulties associated 
with the field theory, it seems fruitful to search for al- 
ternative approaches to the problem which might shed 
some further light on the interesting properties of the 
surface 0i exponents in d = 1. One such alternative is 
provided by the theory of quantum spin Hamiltonians, 
to which we turn in the next section. 



E. Exact Results 

In this section we will derive some exact results for the 
surface 0i exponents in 1 + 1 dimensions at the O and Sp 
transitions. The methods are a straightforward extension 
of the work in 1 28 3^] . The starting point is the following 
set of rules for BARW with m = 2 in 1 + 1 dimensions: 



%A <-> A% with rate D/2 

AA -> 00 with rate A (42) 

0A0 <-> AAA and §AA <-> AA% with rate a/2. 



The 1 + 1 dimensional regime which should more 
amenable to field theoretic analysis is when both the 
bulk and surface branching processes are unimportant, 
and hence we should be able to use results derived solely 
from the A + A — > reaction. This is the case for the 
region <r m < a m , critical in d = 1 (see Figure |a). The Sp* 
"transition" has been fully analyzed in Ref. j3^|, which 
predicts a i -1 / 2 decay both in the bulk and at the wall, 
with a density excess at the wall. We note that the crit- 
ical Sp* state in d = 1 can be characterized in two ways: 
as a decaying density, or as a survival probability. If we 
place two particles close together in the bulk at t = 0, 
then, from simple random walk theory, the probability 
these particles are still alive at time t scales as t^ 1 / 2 . 
Hence, in the bulk, these two ways of characterizing this 
phase scale in the same way. However, if the two particles 
are released next to the wall, then it is easy to show (us- 
ing the method of images) that the survival probability 
now decays as t^ 1 [ p7|j3q ] . Therefore, these two charac- 



Note that these rules are fermionic in character (no more 
than one particle per site is permitted) in contrast to the 
bosonic rules employed in the derivation of the earlier 
field theory. The model described in ( |42] ) can be trans- 
formed into a spin picture by writing the configuration 
of a semi-infinite system as a vector \s±, S2, S3, . . .}, where 
Si = 1/2 if the i-th site is empty, and Sj = —1/2 if that 
site is occupied. Hence the system ket is given by 

|i , (*)) = E p « s «>;*)i^»' ( 43 ) 

and the equation governing the time evolution is 

d t \P(t)) = -H\P(t)), (44) 

where, using a representation in terms of Pauli ma- 
trices, and defining rik = (1 — cr|)/2, Vk = 1 — rik, 
s fc = (°fc ± *°fc)/ 2 > we nave @ 
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1 



2 ^ 

k=l 



(D[n k v k+1 + v k n k+ i - s 



k °k+l 



'fc "fc+U 



+2X[n k n k+ i 
= DH SEP + 



44+J) + % E( : - Vk-iv'k+iW 



XH 



RSA 



fc=2 

aff BARW 



(45) 



^E^ EP 

fe=i 



aE^ sa +«E^ arw - 



fc=2 



Here we have used some of the notation of [p8[ , where 
SEP (symmetric exclusion process) refers to the diffusion 
piece, RSA (random-sequential adsorption) to the anni- 
hilation piece and BARW to the branching piece of the 
"quantum Hamiltonian" . Notice that the boundary has 



been included in ([45), since particles may not hop to the 
left of site 1, and the annihilation/branching processes 
have also been restricted to sites 1,2,3.... Hence, the 
above operator Tt governs the evolution of a 1 + 1 dimen- 
sional BARW system without an A — > reaction at the 
boundary. Averages are calculated using the projection 
state (| = E{ Si }(MI> i-e- (?) = <|^TO>- Following 
p8[ , we now define an operator D where 



V 



7-i7o7i72 ■ 



(46) 



with 



72fe-l 



5 [(i+^K-(i 



■i)], 



TS»fc = 5 [(l + iK^i-(l-0]- (47) 

Defining a new "quantum Hamiltonian" as Ji = 
[D~ 1 HV] T , we find 



H = [D-X}T h* ARW 



fe=i 



A] 



E^ EP 

k—1 



[mn 



5-i 5 n 



mv 



>E^ SA 

]. (48) 



where we have used the commutation rules described in 
detail in p8|. Hence, when D = X + a, we have the 
following processes occurring: 

%iA i+1 <-> A$ i+ i rate (A + a)/2, 

* = 1,2,3,..., 
1 = 1,2,3,..., 
i = 1,2,3,..., 
1 = 1,2,3,..., 



AiA i+ i — > 0j0i+i rate A, 
-iA$ i+ i <-> A i - 1 A l A l+1 rate a/2, 
<-> Ai-iAi$ i+ i rate a/2, 
A Ai -> o 0i rate A/2, 

%Ai -> A o 0i rate A/2. 



(49) 



Excepting the boundary terms, we see that the Hamilto- 
nian has been mapped back onto itself. Furthermore, at 
the edge, the particles may only hop from site 1 to site 
but never the other way round. This means that we can 



forget about the 0-th site in exchange for allowing the 
processes A\A 2 «-» Ai0 2 (with rate a/2), and A\ — > 0i 
(with rate A/2). Hence, we see that the new Hamilto- 
nian TL corresponds to the case where [i s ^ 0, with the 
DP processes A «-> A + A and A — > generated on the 
boundary. 

If we choose the initial condition to be an uncorrelated 
state with density 1/2, denoted by |l/2), then the density 
at site k, p k {t), is given by 



Pk{t) = (Kexp(-?#)|l/2}. 



(50) 



Following exactly the procedure in |2^,|39| (starting with 
insertions of the identity operator X>X> - Mnto the RHS of 
( [sb] ) ) , one can straightforwardly show that 



Pk(t) = g[l 



\exp(-nt)\k - l,k)}, 



(51) 



where (0| is the vacuum state (with no particles), and 
fc — 1, fc) is the initial state with only two particles situ- 
ated at sites fc — 1 and fc. However the relation ( |5l| ) is just 
what we wanted to prove: the LHS is the density at the 
fc-th site, whereas the RHS is one-half times the proba- 
bility that a cluster initiated at t — by two particles at 
sites fc — 1 and fc has not yet died out by time t. Accord- 
ing to our earlier analysis, for A < 0, the LHS should 

scale as |A|^ dcnB (far from the wall) or |A|^ 1 - d « M ' (close to 
the wall), and the RHS as |A|' 3seGd (far from the wall) or 
|A|^.«*»d (close to the wall). Thus, at the line D = X + a, 
we have shown the desired result /3p scod = @ x ^ ens (and, 
of course, the bulk result /3 see d = /?dcns)- We note that 
the bulk result was proven in |p8| . and a very similar re- 
sult for A + A — > was derived in |3^] (connecting the O* 
and Sp* "transitions"). Using universality, we postulate 
that the equality between the two surface exponents is 
valid everywhere close to the transition line, and not just 
where D — X + a. 

It is now straightforward to derive the relation 
^i,seed = ^°dens (again at the line D = X + a). One 
simply starts off with the quantum Hamiltonian Ti. and 
then follows the same steps as above. H can then be 
mapped back onto the starting Hamiltonian 7i, meaning 
that the transformation is actually a duality transforma- 
tion. A relation like that in ( |5l| ) can then be derived, 



giving pf p seed - Mi, dons - 

In summary, at the particular line in parameter space 
D = A + a, we have derived some exponent equalities 
which are in full agreement with the simulations to be 
presented in the next section. In particular, we see that 
we have mapped BARW at the special transition onto 
BARW at the ordinary transition in 1 + 1 dimensions 
(and vice- versa), a rather non-trivial procedure. This 
has allowed us to derive some results about the (3\ expo- 
nents (something which seems to be beyond the ability 
of the field theoretic methods at present). Unfortunately, 
as is always the case with exact calculations, the result is 
only derived for one line in parameter space and we have 
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to rely on universality in order to claim that it is valid 
elsewhere close to the transition line. 



V. DPn AND BARW MODELS 

We will now briefly present the specific models and 
boundary conditions used in our numerical simulations. 
We begin with DP, include its generalization to DP2, and 
then comment on how we implement BARW. In all cases 
we include the specific boundary conditions a nd identify 
them according to the classification in sections III and IV . 

For d = 1, bond DP as well as site DP (for which 
the sites percolate instead of the bonds), are contained 
in the Domany-Kinzel model |i(],[4l[]. Each site can ei- 
ther be active or inactive and the probability for site i 
to be updated to state s^t+i at time t + 1 is given by 
P(si,t+i\si-i,t, Sj+i,t)- See Figure ^| for a typical lattice 
configuration and Figure fj] for the update rules. 



space 



time 
t 

t+1 



o o; 
p.; 

A • 

D 



A 

A 

A 



-Qh 



FIG. 6. Directed Percolation in terms of the Do- 
many-Kinzel model, where time flows vertically downwards. 
Black sites are active (A) and white ones inactive (I). The 
state of each site at time t+1 depends on the states of the 
neighboring sites at time t. 



O. 



O 



Q. 



• • 

1 p q 

FIG. 7. Update probabilities for DP in terms of the pa- 
rameters < p, q < 1, where we have q = p(2 — p) for bond 
DP and q = p for site DP, respectively. Probabilities for the 
other configurations follow from left-right symmetry and from 
P(A\...)+P(I\. ..) = !. 

The DP2 model has two symmetric absorbing states in 
which the system can be trapped. It is a special case of a 
generalized Domany-Kinzel model (DPn) introduced by 
Hinrichsen [ [l2"[ , where each site can be either active or in 
one of n inactive states. For n = 1 the update rules are 



identical to those of the Domany-Kinzel model in Fig- 
ure 0, but for n > 2, the distinction between regions of 
different inactive states is preserved by demanding that 
they are separated by active ones. An example of a DP2 
cluster is shown in Figure |^b, where we have also shown 
an ordinary DP cluster in Figure 0a for comparison. In 



1 + 1 dimensions, DP2 belongs to the BARW universalit 



class, and the update probabilities are given in Figure 



(b) 



I 




FIG. 8. (a) A DP cluster and (b) a DP2 cluster both grown 
from a single seed in the bulk. 
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FIG. 9. Update probabilities for DP2: black sites are 
active (A), whereas white and grey sites are in the inac- 
tive states Ii and I2, respectively. Probabilities for the 
other configurations follow from left-right symmetry and from 
P(A I ...)+P(Ji I ...) + P(h I ...) = 1. 

The easiest way of introducing a boundary into DP 
and DP2 is simply to cut off the lattice. This is equiv- 
alent to introducing boundary sites which are forced to 
be in one of the inactive states. We will refer to this 
case as the inactive boundary condition (IBC) and we 
choose inactivity of type 1 to the left of the boundary, 
see Figure nq. Apart from imposing the state of these 
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sites within the wall, the sites at the wall and those in 
the bulk are updated by the rules in Figures @ and |. 




FIG. 10. DP2 with an inactive boundary condition (IBC), 
corresponding to the special (Sp) universality class. 

Next we consider the reflecting boundary condition 
(RBC) where the wall acts like a mirror so that the sites 
within the wall are always a mirror image of those next 
to the wall, see Figure [fl]. For DP2, one can see that 
there is a qualitative difference between the IBC and the 
RBC. For the latter, regions of type-2 inactivity can get 
trapped at the wall and the only way for these regions to 
disappear is to wait for the cluster to return, whereas for 
the IBC such regions are never trapped. 




FIG. 11. DP2 with a reflecting boundary condition (RBC), 
corresponding to the ordinary (O) universality class. 

We now consider the active boundary condition (ABC) 
where the sites within the wall are forced to be active, 
see Figure [l2| In this case the cluster will never die com- 
pletely as the wall will always be active and can always 
induce new clusters. Nevertheless, by defining the sur- 
vival time of a cluster as the point in time when the 
system has no activity apart from within the wall itself, 
we can define the same exponents for the ABC as for 
the other boundary conditions. We have, however, not 
studied this boundary condition in any detail but merely 
mention it here for completeness. 

We can now discuss the relation between the above 
boundary conditions and our previous classification of 
the universality classes at the boundary for BARW in 



1 + 1 dimensions. The key feature is whether the sym- 
metry between the two absorbing states in the bulk is 
preserved at the surface. In terms of the DP2 model, the 
IBC model respects this symmetry and hence it belongs 
to the special (Sp) universality class, whereas the RBC 
model does not respect this symmetry, and hence belongs 
to the ordinary (O) universality class. Furthermore, the 
ABC model clearly belongs to the extraordinary (E) uni- 
versality class. Hence we see that by using the IBC, RBC, 
ABC classification all the previously discussed boundary 
BARW transitions in 1 + 1 dimensions can be accessed. 




FIG. 12. DP2 with an active boundary condition (ABC), 
corresponding to the extraordinary (E) universality class. 

Furthermore, let us note that for DP the classification 
of the IBC and RBC is somewhat different. In 1 + 1 di- 
mensions DP probably does not support the special tran- 
sition, and since no symmetry is broken by the RBC, both 
the IBC and RBC will belong to the ordinary transition 
universality class ]li| , p2| . 

We have also performed simulations for a lattice 
BARW model with IBC and RBC boundary conditions. 
For BARW we initially placed two particles at the two 
sites closest to the wall. The one-dimensional BARW 
model is then implemented with "dynamic branching" , 
which means that the branching of one particle into three 
particles occurs randomly to either the left or the right 
of the particle |^] . The BARW model is expected to be 
in the same universality class as DP2 and our results are 
in agreement with this for both bulk and surface quan- 
tities (using the IBC and RBC). Hence, in the following 
section, we will only discuss the results for DP2 since 
this is also the model to which we devoted most of our 
simulations. 



VI. NUMERICAL RESULTS 

We have studied DP2 in detail using Monte Carlo sim- 
ulations in 1 + 1 dimensions. The wall is placed at i = 
and we use an initial configuration with one active site 
at i = 0, with the sites i > being in the inactive state 
h ■ Thus the absorbing state corresponds to the situation 
where all sites are in the inactive state I\. The system 
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is evolved according to the DP2 rules (see Figure^), and 
we typically average over 10 5 independent clusters in or- 
der to reduce the error bars to a few percent. Using the 
notation of Figure ^, we have carried out simulations for 
q = p at the critical probability p c , where we have used 
the estimate p c — 0.5673 fl~2]| . 

In these simulations, starting from a seed on the wall, 
we measure the survival probability P\ (t) , the activity in 
the bulk Ni(t) and at the wall Ni t \(t), the average spread 
of the cluster (a; 2 (i)), and the probability p±(s) to have a 
cluster of size (mass) s, all at criticality (T^|2^]. Further- 
more, by averaging over surviving clusters only (denoted 
by an over-line), we measure the surviving bulk activ- 
ity Ni(t) and the surviving wall activity N\ t i(t), again 
starting from a seed on the wall. 

First we performed simulations for DP2 without a wall 
and obtained results for the exponents in complete agree- 
ment with those in jl^] . Our results are listed in Tables || 
andf| There are several estimates available for the bulk 
exponent /3dcns (= Aseed) p3| ]- In the following we will 
use the estimate f3 dcns = 0.922(5) @. 

We now list some exponent relations used to extract 
the exponents from our numerical simulations ]45| . All 
the relations given below are valid for both the IBC (spe- 
cial) and RBC (ordinary) DP2 transitions, and hence 
these labels will be suppressed from now on. The prob- 
ability for a cluster grown from a seed on the wall still 
to be alive at time t is given by Eq. (^4|). At criticality 
(A = 0) it has the following behavior 

Pi(t) ~ t- 5l '»°° d , (52) 

with the exponent 

5l,seed = Pl,aeed/v\\- (53) 

Hence, the probability of growing a cluster which lives 
exactly t time steps behaves as Px(t) ~ t~ 1 ~ Sl - Beed . Away 
from criticality it is straightforward to obtain the average 
cluster lifetime of finite clusters from (|2~i|). One obtains 

<t)~|Ar\ (54) 

with the exponent 

T\ = f || - ^l,seed- (55) 

The average number of active sites at criticality, av- 
eraged over all clusters, is obtained by integrating the 
density (|25|) over space, and one arrives at 

Nx(t)~t K \ (56) 

with 

K l = dX - Helens - £l,seed) (57) 

where we have introduced the cluster envelope or "rough- 
ness" exponent x = v ^l v \\ (= }_L z )i ano - the notation 
<5dcns = Pdcns/v\\. Note that ( p7| ) corresponds to the 



hyperscaling relation (p0[), a fact which follows from 
(si(i)) = f dt' Ni(t'), and the relation 71 = i/m(1 + Ki). 
By integrating the density on the wall (p6|), we obtain 
the average number of active sites at criticality on the 
wall 

JV u (t)~f ,1 .\ (58) 

with 

=(d- 1)X - <5l,dcns - <5l,sccd, (59) 

and where <5i.dcns = /8i,densA / ||- Note also that ( |59| ) cor- 
responds to the hyperscaling relation ( |3^ ) at criticality, 
since (si,i(t)) = J dt' N 1A (t r ), and 71,1 = v\\ (1 + 

Alternatively, by averaging only over clusters which 
survive to infinity (denoted by an over-line), we obtain 

Nl(t) ~ t~, (60) 

where 

Ttl = dX - ^dcns- (61) 

The activity on the wall for surviving clusters reads 

iVi7(i) - 1~, (62) 

with the exponent 

WJ = (d - l)x - <5l,dens- (63) 

Simulations in 1 + 1 dimensions thus directly yield <5i,dcns- 
The average position of activity follows from ( p5| ) 

(x 2 )~t 2x , (64) 

where x is the distance from the seed and the average is 
taken over all active points at a given time. 

For further confirmation of our numerical data we also 
considered the cluster size distributions at criticality. In 
the bulk the typical cluster size s of finite clusters scales 
as volume times density, i.e. 

a ~ £l£||n(A) ~ |A|-V-, (65) 

with 

1/(7 = dvj_ + Z/|| - /?dons- (66) 

From the lifetime survival distribution (|52]), it is then 
straightforward to obtain the probability Pi(s) to have a 
cluster of size larger than s, for clusters started from a 
seed on the wall. Using the fact that the lifetime is set by 
the parallel correlation length, t ~ £m ~ |A| _i/ ii, we see 
that the typical cluster size and lifetime are connected 
by 

s^t 1 /^". (67) 

Hence we obtain P x (s) - P x (i - s u \\ a ) - S -Pi,^°. 
Thus, we eventually obtain the probability p\{s) to have 
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a cluster of exactly size s, pi(s) 
result 



where 



Mi 



1 



Pl (s)~ s -^, 

Pl.seed 



-dP\(s)/ ds, with the 



(68) 



dv± + v\\ — (3d 



(69) 



Similarly, the cluster size distribution on the wall due to 
a seed located at the wall can also be obtained. In this 
case the typical cluster size of finite clusters is 



where 



Swall-^VltA^IAr 1 /- 1 , 



l/(7i = {d— l)vx + V\\ - /3l,dens- 



The resulting distribution reads 



Pl,l( s wall) ~ S 



Mi,i 
wall " 



with 



A*l,l 



l.sccd 



(d - l)v±_ + - /3i,dens 



(70) 



(71) 



(72) 



(73) 



Note also that many of the scaling expressions given 
above only apply exactly at bulk criticality. Away 
from that point one must also include a scaling func- 
tion. For example, Eq. ( p2| ) is replaced by P\{t, A) = 
lA^i.-cd F{t/\A\- V u), and Eq. (f|) is replaced by 
Pi(«,A) = |A|^ I1 / CT G(s/|A|- 1 /^). 

In Tables | and || we list our estimates for the critical 
exponents for DP2. Our results are in complete accor- 
dance with our theoretical analysis: bulk exponents are 
unaltered whereas the wall introduces two separate sur- 
face exponents. We have also carried out bulk and surface 
simulations for A < and confirmed that our data can 
be collapsed according to an appropriate survival proba- 
bility scaling function [see ( pif ) for the surface case] , using 
our exponent estimates. This numerically confirms the 
validity of the relation 5 = (3/v\\ for the bulk, as well 
as the analogous relations for the both sets of surface 
exponents p6[ . 

By using the explicit definitions of the IBC and RBC, 
we can deduce some further properties of the /3i )See d and 
/3i.dcns exponents. There will be more activity next to 
the wall for the IBC than for the RBC, since the latter 
can have regions of I2 located at the wall. Once created, 
these I2 regions will survive until the activity returns to 
the wall. Thus it follows that /3} B d ^ ns < (3f^ ns . On the 
other hand, the existence of these I2 regions implies that 
the survival probability (^i|) for the RBC will be greater 
than for the IBC, leading to /3{ B ° ed > P?®£ d . Note that 
both our simulations and our previous exact calculations 
show that /3i,seed + /3i,dens is the same for both the RBC 
and IBC. Using a hyperscaling relation [like that in (p2m, 



this implies that 7i^ c = J^ 30 (although both exponents 
defined in this way are negative). We have also studied 
several other boundary conditions and found that these 
give the same scaling behavior as either the RBC or IBC 
depending on whether the above-mentioned I2 regions 
can disappear only at the wall or also in the bulk. 

We can obtain an interesting exponent relation for the 
RBC transition by assuming that the survival probabil- 
ity is dominated by the return to the wall of the cluster- 
envelope, which leads to [15[ 



RBC 

iced 



1 



(74) 



in agreement with our simulation results for the RBC. 
Qualitatively, this means that the I2 regions located at 
the wall determine the scaling since they can only disap- 
pear when the activity returns to the wall. Note that a 
relation of this kind clearly fails for the IBC transition. 
Furthermore, if the cluster lifetime is defined to be the 
return time of the cluster-envelope (i.e. the return time 
of the rightmost active site) to the initial point, then 
we expect clusters defined in this way to have a lifetime 
distribution exponent <5i iSOC d given by Eq. (f74|). This pre- 
diction is in agreement with the simulations in [^ , where 
various models in the DP and BARW classes were stud- 
ied with cluster lifetimes defined in the way described 
above. 

For DP it has been customary for some time to in- 
vestigate whether the critical exponents can be fitted by 
simple rational numbers |^7j . Such a fitting has also been 
tried for bulk BARW with the following guesses in 1 + 1 
dimensions: n = \ — 2<5 = and \ = 4/7 M. These 
estimates lead immediately to 5 — 2/7 (and /3/Vi_ = 1/2, 
7 = vu). It is intriguing to note that our numerical re- 
sults for DP2 in addition suggest that [i\ — 3/2 for the 
IBC and 4/3 for the RBC. From Eq. (f|), it then fol- 
lows that <5i, SO cd = 9/14 for the IBC and 3/7 for the 
RBC. We would need one more relation in order to ob- 
tain the last independent exponent, which we can take 
to be v\\ . In fact, we observe numerically that the rela- 
tion 2i^|| — /3i iS ood — /3i,dcns = 3 is valid to within one 
percent. From these observations the remaining DP2 
exponents follow: /3dcns = /iced = 12/13, v\\ = 42/13 
and = 24/13. Furthermore, /3i jSC cd = 27/13 and 
/3i,dens = 18/13 for the IBC, and vice versa for the RBC. 
However, at present we have no understanding of these 
possible exact values for the 1 + 1 dimensional exponents. 
BARW is certainly not conformally invariant and conse- 
quently until some theoretical framework is proposed to 
explain why these exponents could be rational numbers, 
numerical coincidence remains a distinct possibility. 



VII. CONCLUSION 

In this paper we have presented a study of critical sur- 
face effects in systems with non-equilibrium phase tran- 
sitions. In particular we have focused on the DP and 
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BARW universality classes, where we have put forward 
a unified presentation involving mean-field, scaling, field- 
theoretic and exact methods. Furthermore, many of our 
theoretical conclusions have been backed up by large- 
scale Monte-Carlo simulations. 

Nevertheless, there are still a number of open ques- 
tions. In particular our understanding of surface BARW 
in 1 + 1 dimensions is hampered by the fundamental prob- 
lems of the field theory, which mean that the boundary 
(and bulk) transitions occurring at a m = <7 m . critical re- 
main difficult to treat using RG methods. Furthermore, 
we have concentrated on the ordinary and special transi- 
tions in 1 + 1 dimensional BARW — there will most likely 
also be interesting behavior for the extraordinary transi- 
tion. 

Finally, our most important result is the existence of 
two independent surface exponents: /3i,dens and /3i, S eed 
for surface BARW (and DP2). This certainly distin- 
guishes DP from BARW, since for the former case, 
/3i,seed = /?i,dcns- For the 1 + 1 dimensional BARW case 



on the other hand, we have used exact techniques to link 
the surface exponents at the ordinary and special tran- 
sitions, giving /3° see <j = (3f p dons and /?f p spcd = ft° dcns . It 
would certainly be instructive to rederive these results 
from a field theoretic perspective, but this is beyond the 
scope of the present paper. 
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TABLE I. Critical exponents obtained from our DP2 simulations. For comparison we also list the exponents for DP in the 
bulk and with an IBC wall The <5i lSOC d exponent is obtained from Eqs. (^) and ([57]). Exponents without the '1' 

subscript refer to the bulk. 





DP 


DP (IBC) 


DP2 


DP2 (IBC) 


DP2 (RBC) 


Kl 


0.473 14(3) 


0.473 14(3) 


0.288(5) 


0.287(2) 


0.285(2) 


<^dcns 


0.159 47(3) 


0.159 47(3) 


0.287(5) 


0.288(2) 


0.291(4) 


Pd ens 


0.276 49(4) 


0.276 49(4) 


0.922(5) 


0.93(1) 


0.94(2) 


K. 


0.313 68(4) 




0.000(2) 






Kl 




0.0496(3) 




-0.354(2) 


-0.141(2) 


$l,seed 




0.4235(3) 




0.641(2) 


0.426(3) 


Pi, seed 




0.7338(1) 




2.06(2) 


1.37(2) 


$l,dens 




0.4235(3) 




0.415(3) 


0.635(2) 


/?l,dens 




0.7338(1) 




1.34(2) 


2.04(2) 


2X 


1.265 23(2) 




1.150(5) 


1.150(3) 


1.152(3) 


X 


0.632 61(2) 




0.575(3) 






VII 


1.733 83(3) 




3.22(3) 






v± 


1.096 84(2) 




1.84(2) 
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TABLE II. Critical exponents for cluster lifetime (^) and mass distributions (|§), (||), @, and @. For comparison we 
also list the exponents for DP in the bulk and with an IBC wall |ilj||l8||l7| . We also give the exponents for the average lifetime 
(p4j), and average cluster sizes (|29|) and (blh , obtained from the scaling relations. 





DP 


DP (IBC) 


DP2 


DP2 (IBC) 


DP2 (RBC) 


^seed 


0.159 47(3) 




0.290(5) 






^l,seed 




0.4235(3) 




0.646(3) 


0.425(3) 




1.108 25(2) 




1.225(5) 






Mi 




1.2875(2) 




1.500(3) 


1.336(3) 


Mi,i 


1.189 72(6) 


1.7337(2) 


1.408(5) 


2.05(5) 


2.15(5) 


T 


1.457 34(7) 




2.30(3) 






Tl 




1.0002(3) 




1.16(4) 


1.85(4) 


7 


2.277 69(4) 




3.22(5) 






7i 




1.8207(4) 




2.08(4) 


2.77(4) 


71,1 


1.180 85(4) 


0.2664(3) 


1.38(3) 


(<0) 


«0) 
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